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PLANE-STRESS  ANALYSIS  OF  AN  EDGE-STIFFENED 
RECTANGULAR  PLATE  SUBJECTED  TO  I/DADS 
AND  TEMPERATURE  GRADIENTS 


By  Charles  Iobove,  Dalpat  Ranchal,  and  Frank  Dunn 

SUMMARY 


(a plane-stress  analysis  is  presented  for  an  isotropic  or  orthotropic 
elastic  rectangular  plate  hounded  hy  four  edge  stiffeners  and  subjected  to 
prescribed  temperature  distributions  and  boundary  loads.  The  stiffeners 
are  assumed  to  have  zero  bending  stiffness  but  finite  extensional  stiff¬ 
ness.  They  may  be  either  uniform  or  tapered  in  such  a  way  as  to  develop 
a  constant  (or  non-constant  but  prescribed)  cross-sectional  stress.  In 
the  latter  case  the  required  form  of  taper  is  determined  along  with  the 
stress . 

The  solution  is  "by  means  of  Fourier  series.  Its  convergence  and 
feasibility  are  tested  by  means  of  two^ numerical  examples ,  one  a  thermal- 
stress  problem,  the  other  a  " shear- lag  1  problem^  ^  ^ 

INTRODUCTION 


A  plane-stress  analysis,  by  means  of  Fourier  series,  is  presented  for 
an  elastic  rectangular  plate  bounded  by  four  edge  stiffeners  and  subjected 
to  prescribed  loads  and  temperature  distributions .  The  plate  may  be  iso¬ 
tropic  or  orthotropic,  with  elastic  constants  that  are  independent  of 
position  and,  if  orthotropic,  with  axes  of  elastic  symmetry  parallel  to 
the  edges.  The  four  edge  stiffeners  are  assumed  to  have  zero  bending  stiff¬ 
ness  but  finite  extensional  stiffness  and  to  be  integrally  attached  to  the 
plate  along  their  originally  straight  centroidal  axes.  The  stiffeners  are 
either  uniform  or  tapered  in  such  a  way  as  to  have  constant  (or  non-constant 


SERIES  EXPANSIONS  FOR  PRESCRIBED  LOADS  AND  THERMAL  STRAINS 


„  ^results  °f  th®  present  analysis,  to  be  discussed  shortly,  consist 
omiulas  for  the  stiffener  and  plate  stresses  in  terms  of  the  given  load- 

^oSno+  erml"s^ra^  distributlon*  However,  the  loading  and  thermal  strains 
do  not  appear  explicitly  in  these  formulas;  it  is  rather  the  Fourier  coef¬ 
ficients  of  these  quantities  that  are  required.  In  anticipation  of  this  re- 
quirement,  it  is  assumed  that  the  given  distributed  loadings  can  be  expressed 
in  Fourier  series  of  the  following  form,  with  known  coefficients: 


,  .  N 

Nqty)  =  n21  B^  sin  (mty/b) 

~  n=l  Bn  sin  (n*yA>) 

M 

N3^X^  =  mil  Bmsin  (mitx/a) 

M 

Vx>  =  mil  Bm'  sin  W) 

,  ,  N 

ql'y'  ~  nio  %  cos  (nny/b) 

(y )  =  n|Q  <£  cos  (njy/b) 

,  \  M 

q3X  ~  mio  cos  (nutx/a) 

M 

=  mi&  ST  cos  (mnx/ a ) 


(3) 


S+^iiar'ly *  there  are  any  discontinuities  in  thermal  strain  between  the 

S^1f;en!r®  and  the  plate  edges,  these  are  assumed  to  be  known  in  the  form 
oi  the  following  Fourier  series : 
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el(y)  -  ey(0,y)  =  Jx 

N 

e2(y)  -  ey(a,y)  =  ^ 

,  ,  M 

e  (x)  -  e  (x,0)  =  S 

5V  '  x  m=l 

M 

e,  (x)*  -  e  (x.b)  =  Z 

x'  ’  m=l 


TT  sin 
n 

(nity/b) 

for 

T"  sin 
n 

(mty/b) 

for 

Tm  sin 
m 

(mitx/a) 

for 

Tn"  sin 
m 

(rmtx/a) 

for 

0  <  y  <  b 

0  <  y  <  Id 

(5) 

0  <  x  <  a 

0  <  x  <  a 


Finally,  d2e  /dx2  +  52e  /cty-2  is  assumed  to  be  representable  by  the 

y  x 

following  series  in  the  open  region  0  <  x  <  a,  0<y<b: 


52ey  S2ex 

ax2  ^2_ 


M  H 
m=l  n=1 


imrx 

sin  -  sin 

a 


n  Tty 
b 


(6) 


Finite  upper  limits  M  and  N  are  shown  for  the  summation  indexes  in 
equations  (3)  to  (6)  in  expectation  of  the  fact  that  it  will  normally 
be  necessary  to  use  truncated  rather  than  infinite  series  for  practical 
computational  reasons. 


In  the  analysis  of  the  constant -area-stiffeners  case  in  appendix  B 
the  Fourier  coefficients  of  the  loading  and  thermal  strains  will  appear  in 
certain  groups.  These  groups  are  designated  by  the  letters  R^  R^11  , 

R,,n  for  convenience  and  are  defined  by  equations  (B 66)  and  (b67)* 
m 

The  Fourier  coefficients  appearing  in  equations  (3)  to  (6)  and  required 
for  the  evaluation  of  the  Rr .  R",  R,M.  R,,n  can  be  determined  from  the  usual 
definitions.  For  example. 
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where  §,  Q  is  Kronecker's  delta.  Integration  by  parts  in  the  last  equation 

gives  the  following  alternate  formula  which  permits  T  to  be  evaluated 

inn 

from  the  first  derivatives  of  e  and  e  instead  of  the  second  derivatives: 

y  x 


b  a 


mrc  4 

a  ab 


y  mjtx  .  nicy  ^ 

nr cos  ~ sin  b  dy 


0  0 


a  b  -n 

nn  k  T  P  <*x  nny  . 

'“Sb  J  J  cos  b  S1 
0  0 


iny  .  mrtx  _ 

rf-  sin  -  dy  dx 

b  a 


Equation  (9)  may  be  used  for  discontinuous  e^  or  e^  provided  that  Sey/ck 

and  (5e^/by  are  regarded  to  be  infinite,  in  the  manner  of  the  Dirac  delta 

function^  at  the  loci  of  points  of  discontinuity.  If  e_  and  e__  are  con- 

tinuous  in  the  region  0  <  x  <  a,  0<y<b,  further  integration  by  parts 
gives 


?mn  =  "  IT  SF  f  [ey(a>y)  cos  -  ey(°>y)l  sin  SSL  dy 
0 


/HITC  \  2  4 

{  a  ;  ab 


b  a 


0  0 


ey(x*)y)  sin  sin  £x  cly 


nir  4 


b  ab 


7C  /  [eY(x^)  cos  -  e  (x,0)]  sin 


a  b 


-  (— )2  — 

'  b  '  ab 


ex(x,y)  sin  SSL  sin  SSL  a.y  a* 


0  0 


RESULTS  FOR  CONSTANT-AREA  STIFFENERS 

Series  for  the  plate  and  stiffener  stresses.-  The  analysis  in  appendix 

B  gives  the  plate  stress  resultants  and  stiffener  forces  in  the  form  of 

series.  Equations  (B2l)  to  (B27)  give  the  plate  stress  resultants  N  ,  N  , 

y  xp 

and  Nxy,  equations  (BI7)  the  stiffener  forces.  It  should  be  noted  that  the 

series  given  for  N  in  the  interior  of  the  plate,  equation  (B2l),  is  not 

»/ 
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valid  along  the  edges;  separate  series,  equations  (B22)  and  (B23),  are 
given  for  evaluating  N  along  the  edges  x  =  0  and  x  =  a;  along  the . other 
two  edges,  y  =  0  and  y^  =  h,  is  equal  to  N^(x)  and  N^(x)  respectively. 

Analogous  remarks  apply  to  N  .  Similarl y  the  series  for  the  stiffener 
forces,  equations  (B17)>  are  not  valid  at  the  ends  of  the  stiffeners; 
this  imposes  no  great  hardship,  however,  since  the  stiffener  end  forces 
are  known  as  part  of  the  given  loading . 


Evaluation  of  series  coefficients.  -  In  order  to  use  these  series  for 
numerical  calculation  of  stresses,  one  must  first  evaluate  the  coefficients 
c'  c"  g'  g",  c  ,  g  ,  j  .s’  s",  s'",  and  s""  appearing  in  them.  The 

cn’  n*  &m’  mu’  Bmn’  n*  n*  m  m 

first  four  groups  of  coefficients,  namely  e^,  c^,  g^,  g^,  are  the  key  to  all 
the  others,  so  their  evaluation  will  he  discussed  first. 

The  c' .  c".  g' ,  and  g"  are  defined  by  the  system  of  equations  (B6l)  to 
n’  n’  V  m  .  . 

(b64)  and  can  be  determined  by  solving  these  2N  +  2M  equations  simultaneously 

for  the  c*  c".  a '  g"  and  noting  the  definitions  in  equations  (B65).  As  an 
n’  n’  &nr  °m  „  ,  .  „ 

alternative,  equations  (B70)  may  be  solved  simultaneously  for  the  gffl  and  g^, 

the  c'  and  c"  are  then  obtained  directly  from  equations  (B68) .  This  alterna¬ 
tive  is  preferable  because  it  requires  the  solution  of  only  2M  simultaneous 
equations,  regardless  of  how  large  a  value  is  selected  for  N. 

If  the  structure,  loading  and  thermal  strains  are  symmetrical  about  the 
centerline  y  =  b/2,  considerable  simplification  results.  The  c^,  cn,  g^, 

and  e"  are  then  defined  by  equations  (B78)  and  the  simultaneous  system  of 
M  +  N“+  1  equations  (B6l* )  to  (b63*).  As  a  preferable  alternative  to  equa¬ 
tions  (b61')  to  (B63' ),  the  M  equations  (B70')  may  be  solved  simultaneously 
for  the  g^,  after  which  the  g^,  c^,  and  c^  are  obtained  directly  from  equa¬ 
tions  (B78)  and  (B68) .  In  this  alternative  the  size  of  N  again  does  not 
influence  the  number  of  equations  that  have  to  be  solved  simultaneously. 

If  both  y  =  b/2  and  x  =  a/2  are  axes  of  symmetry  of  the  structure, 
loading,  and  thermal  strains,  then  equations  (B79),  (B6l  ),  (B63  J  may 
be  used  to  obtain  the  c^,  c",  g^,  g^,  where  (B6l  )  and  (B63  )  represent 

(M  +  N  +  2) /2  simultaneous  equations.  The  quantities  may  also  be  deter¬ 
mined  from  fewer  (namely  (M  +  l)/2)  simultaneous  equations  by  using  the 
system  (B70m)  to  solve  for  the  g^  and  then  obtaining  the  remaining  coef¬ 
ficients  directly  from  equations  (B79)  and  (B6l" ) . 

With  the  ef,  c",  g',  g”  known,  equations  (B59)>  (b60),  (B58),  and 
n  n  m  m 

(B56)  will  furnish  the  remaining  coefficients  directly. 
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ratios  with  each  other  ,  the  solution  takes  on  a  much  simplified  form 
characterized  by  the  fact  that  it  is  no  longer  necessary  to  solve  simul¬ 
taneous  equations.  The  quantities  c',  c",  g*  and  g"  are  defined  by 

n’  n’  Dnr  ^m  J 

equations  (b8o)  directly  on  this  limiting  case. 


In  equations  (b8o)  the  coefficients  of  R' /(A..E  )  etc  M  are  in  error 
by  terms  of  the  order  of  (a^-^AjE  J"1.  A  more  accurate  solution,  in 

which  these  coefficients  are  correct  to  terms  of  the  order  of  (a3Ei:LA  E  )-1, 


is  represented  by  equations  (B8l)  to  (b84) .  This  more  accurate  solution 
still  retains  some  of  the  simplicity  of  equations  (b8o),  in  the  sense  that 
it  too  does  not  require  the  solution  of  simultaneous  equations. 


Numerical  examples.  -  The  foregoing  results  were  used  to  obtain 
numerical  stress  data  for  two  illustrative  problems ,  one  a  thermal-stress 
problem  involving  non-uniform  temperature  distribution  without  any  applied 
loads,  the  second  a  "shear-lag"  problem  involving  the  diffusion  of  loads 
from  the  stiffener  ends  into  the  plate.  The  two  problems  are  shown 
schematically  in  figure  4. 

In  both  problems  the  plate  is  square  (b  =  a)  and  isotropic,  with 
Young fs  modulus  E,  Poisson’s  ratio  v,  and  thickness  h.  The  four  stif¬ 
feners  are  assumed  to  be  identical  and  to  have  the  same  Young's  modulus 
as  the  plate.  The  symbol  A  will  designate  the  common  values  of  A  ,  A  , 

Aj,  and  A^ .  12 

In  the  thermal  stress  problem  (figure  4(a)),  the  stiffeners  are 
at  a  uniform  temperature,  say  zero,  while  the  plate  has  a  "pillow-shaped" 
temperature  distribution  of  the  form  0  sin  (jtx/a)  sin  (jty/a);  thus  Q 
denotes  the  temperature  rise  of  the  plate  center  relative  to  the  edges. 

The  symbol  a  will  denote  the  coefficient  of  thermal  expansion  of  the 
materials . 

In  the  shear-lag  problem  (Figure  4(b)),  the  temperature  is  uniform 
and  the  loading  consists  of  identical  tension  loads  of  magnitude  P  applied 
to  the  end  cross  sections  of  the  stiffeners . 

In  these  problems  the  structure,  loading,  and  thermal  strains  are 

symmetrical  about  both  centerlines,  namely  x  and  y ■=  a/2.  Hence  equations 

(B70")  may  be  used  for  their  analysis.  In  these  equations  R'  and  R,n 

n  m 

take  the  following  forms,  according  to  equations  (B66)  and  ‘  (B67), for 
the  thermal-stress  problem: 


More  precisely,  if  the  dimensionless  quantities  (a^^A^^) with 

i  =  1,  2,  3,  4,  are  assumed  to  approach  zero  while  maintaining  constant 
ratios . 
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R'  =  - 
n 


M 

Z_  K 
m=l  mn 


Rl!  r 


m 


'  N 

Z-  K 
n=l  mn 


where 


K  =  mnjt2T  /(a2E  ) 

mn  mn'  mn7 


E  =  (m2  +  n2f  n4/(a4Eh) 
mn 


'2Q07t2/a2  when  m  =  n  =  1 


mn 


0  otherwise 


Thus 


R* 

n 


EhO0/2 


for  n  =  1 
for  n  f  1 


{Eha0  /2  for  m  =  1 
0  for  m  )  1 


For  the  shear- lag  problem,  on  the  other  hand, 


(ID 


R*  =  R,M  =  kv/a.  (12) 

n  m 


for  all  (odd)  n  and  m. 

The  results  of  the  calculations  are  shown  in  the  figures,  starting 
with  figure  5,  through  dimensionless  plots  of  the  plate  stress  resultants 

N  .  N  .  and  N  as  functions  of  x  for  fixed  values  of  y,  and  the  stiffener 
x*  Y  xy 

force  (x ) ,  (in  view  of  the  symmetry  which  exists  about  the  plate  diag¬ 
onals  in  these  examples,  as  well  as  about  the  plate  centerlines,  the  graphs 
of  p  (y),  P-^y),  and  p^(x)  are  a11  identical  with  ttie  SraPh  for  P3VX)#) 

The  results  shown  are  based  on  calculations  with  M  =  N  =  ?9  •  As  a  check 

on  the  convergence,  the  calculations  were  repeated  with  M  -  99  y  N  =  79 ; 

the  results  agreed  with  the  previous  ones  to  the  extent  that  any  differences 

would  be  indiscernible  on  the  graphs  shown.  The  dotted  curves  in  some  of 

the  figures  represent  approximate  solutions  obtained  by  using  equations  (BoO)  and 

(B8l)  to  (B84),  which  apply  when  the  stiffener  cross-sectional  areas  are 
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large  compared  to  the  plate  cross-sectional  area.  It  is  seen  that  the 
large -stiffener  approximation  is  quite  good  for  the  plate  stresses  when 
the  area-ratio  parameter  4ah/(jr2A)  is  0.1  or  smaller;  it  is  of  course 
exact  when  4ah/(jt2A)  =0. 

The  thermal- stress  results  indicate  finite  shear  stress  in  the  corner 
of  the  plate  for  all  values  of  the  parameter  4ah/(rt2A) .  However,  the  shear- 
lag  problem  exhibits  infinite  corner  shear  stress  for  the  same  values  of  the 
parameter  4ah/(jt  A).  In  an  actual  structure,  finite  stiffener  bending  stiff¬ 
ness  (in  particular,  rigid  connections  where  the  stiffeners  meet  at  the  corners) 
would  lead  to  zero  shear  stress  at  the  corner  and  therefore,  in  a  neighborhood 
near  the  corner,  to  shear  stresses  which  might  be  considerably  different  from 
those  which  the  present  analysis  predicts. 


Accepting  the  premise  of  zero  bending  stiffness  of  the  stiffeners  or 
hinged  connections  where  the  stiffeners  meet,  it  is  interesting  to  note  that 
in  the  thermal-stress  problem  the  maximum  value  of  N  in  the  plate  does 

not  always  occur  at  the  corner  (see,  for  example,  figure  5c). 

.  .  ^  case  4ah/(jr2A)  =  0  represents  the  limiting  case  in  which  the  sheet 

is  infinitely  thin  by  comparison  with  the  stiffener  cross-sectional  dimen- 
sions.  In  this  case  the  end  loads  applied  to  the  stiffeners  in  the  shear- 
lag  problem  should  be  transmitted  unchanged  throughout  the  lengths  of  the 
stiffeners.  Thus  in  figure  6(a)  the  value  of  P  (x)/P  should  be  1.0  for  all 

x/a,  and  N^A/ (hP)  should  similarly  be  1.0  for  all  x/a  along  the  line  y/a=0.0. 

The  deviations  from  1.0  shown  in  figure  6(a)  are  due  to  incomplete  convergence 
associated  with  the  use  of  finite  values  of  M  and  N. 

The  calculations  were  performed  on  the  IBM  7O7O  computer  at  the 
Syracuse  University  Computing  Center  and  required  a  total  of  50.5  minutes 
for  all  the  results  shown.  The  simultaneous  equations  (B70")  were  solved 
by  the  iterative  procedure  of  reference  2.  The  plate  stresses  were  com¬ 
puted  at  x/a  and  y/a  intervals  of  0.1  in  order  to  plot  the  curves  shown. 


RESULTS  FOR  CONSTANT-STRESS 
(OR  PRESCRIBED  STRESS)  STIFFENERS 


For  efficient  design,  it  may  be  desired  to  taper  the  stiffeners  so 
as  to  achieve  constant  cross-sectional  stress  along  the  length  of  each  one 
In  appendix  C  a  generalization  of  this  condition  is  considered  in  which  it 
is  assumed  that  the  tapers  are  such  as  to  produce  prescribed,  but  not 
necessarily  constant,  cross-sectional  stress  variations  along  the  lengths 
of  the  stiffeners.  The  symbols  cr-^y),  a£(y),  cr  (x)  and  cr^(x)  represent 

the  prescribed  stresses  (positive  for  tension)  in  the  stiffeners  located 
at  x  =  0,  x  =  a,  y  =  0,  and  y  =  b  respectively. 

In  this  case  it  turns  out  once  again  to  be  unnecessary  to  solve 
simultaneous  equations  in  order  to  compute  the  plate  stresses.  Equations 
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(B33),  (B34),  (B36),  and  (B3T),  used  in  conjunction  with  equations  (C2), 
give  the  c ' ,  c",  g',  and  g"  directly;  equations  (BoO) ,  (B58),  and  (B56) 

then  furnish  the  j^,  and  equations  (B2l),  (B24),  and  (B27)  the  plate 
stresses • 

The  lengthwise  variations  of  stiffener  cross-sectional  area  needed  to 
produce  the  prescribed  stiffener  stresses  can  be  determined  from  equations 
(C3)  in  conjunction  with  (C4)  to  (C7)* 

CONCLUDING  REMARKS 

r — -  /}// 

| A  plane-stress  analysis ,  based  on  Fourier  series,  has  been  presented 
for  the"^tressesTn"a  linearly  elastic  isotropic  or  orthotropic  rectangular 
plate  bounded  by  four  edge  stiffeners  and  subjected  to  prescribed  tempera¬ 
ture  distributions  and  boundary  loadings.  The  stiffeners  are  assumed  to 
be  either  uniform  or  tapered  in  such  a  way  as  to  develop  constant  (or  non¬ 
constant  but  prescribed)  stresses.  The  convergence  and  feasibility. of  the 
analysis  have  been  tested  and  found  to  be  satisfactory  in  two  numerical 
examples  for  the  case  of  uniform  stiffeners. 

The  present  analysis  differs  from  previous  thermal-stress  analyses 
of  rectangular  plates  by  the  incorporation  of  edge  stiffeners.  It  differs 
from  previous  "shear- lag”  analyses  by  the  avoidance  of  extreme  assumptions 
regarding  the  plate  stiffness  or  plate  normal  stresses  in  what  is  usually 
called  the  transverse  direction. 

However  the  present  analysis  is  also  characterized  by  the  assumption 
of  zero  bending  stiffness  for  the  stiffeners  and  the  absence  of  boundary 
conditions  oh  displacement.  The  removal  of  these  two  limitations  would, 
it  is  felt,  be  a  worthwhile  objective  of  further  research.  Regarding  the 
first  limitation  in  particular,  the  inclusion  of  finite  stiffener  bending 
stiffness  and  joint  rigidity  where  the  stiffeners  meet  appears  to  be 
necessary  in  order  to  obtain  realistic  estimates  of  shear  stress  in  the 
neighborhoods  of  the  corners.  I  / 

— u  /  4 


APPENDIX  A 


SYMBOLS 


Remarks:  (i)  The  subscript  1,2,3  or  4  on  a  symbol  for  a  stiffener- 

related  quantity  identifies  the  stiffener  location  as  x  =  0,  x  =  a,  y  =  o, 
or  y  =  h  respectively.  (ii)  The  Fourier  coefficients  of  known  quantities  * 
(loads,  thermal  strains),  and  combinations  of  such  coefficients,  are  gener¬ 
ally  represented  by  capital  letters ,  while  the  Fourier  coefficients  of 
initially  unknown  quantities  (e.g.,  internal  stresses)  are  denoted  by  small 
letters . 


a  plate  dimension  in  x  direction;  see  figure  1 

amn  Fourier  coefficients  in  series  expansion  for  the  stress 

function  F(x,y);  see  equation  (Bl6) 

ak’an’am’am'U  purler  coefficients  in  series  expansions  for  F(0,y),  F(a,y) 
F(x,0),  F(x,b)  respectively;  see  equations  (Bl4) 

A1  (y )  (y )  >A  (x)  ,A^  (x)  stiffener  cross-sectional  areas 

A  common  value  of  the  above  when  all  four  stiffeners  are 

identical  and  uniform 


b 


Bm" 

rr  n  m  *  m 


c 

mn 


ci>c2>cyck 

e 

mn 


plate  dimension  in  y  direction;  see  figure  1 

Fourier  coefficients  in  series  expansions  for  N  ,N  , 
respectively;  see  equations  (3)  and  (7)  1  2 

Fourier  coefficients  in  series  expansion  for  N  (x,y); 
see  equation  (B2l)  V 

Fourier  coefficients  in  series  expansions  for  N  (0,y) 
and  Ny(a,y)  respectively;  see  equation  (B22)  y 

Cn  C2  cn  C2 

plate  compliances  defined  by  equations  (2) 

Fourier  coefficients  in  series  expansion  for  dSp/cbc^; 
see  equation  (B28) 
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e1(y),  e2(y),  e^(x) ,  e^(x) 


Stiffener  thermal  strains;  see 
figure  2 


ex(x,y),  e  (x,y) 


El>  E2’  e5>  e4 


plate  thermal  strains;  see  figure  2 

C 2 (mjt/ a)^+(C^-2Cj)  (mn/ a ) 2 (n Jt/b ) 2+C ± (mt/b ) ^ 

value  of  above  with  m  =  1  and  n  =  1 

Young* s  moduli  of  stiffeners 

Young* s  modulus  for  stiffeners  and  isotropic 
plate  when  all  have  the  same  Young* s  modulus 


F(x,y) 


stress  function  for  plate;  see  equation  (b4) 


r<  ff" 

’m' 


Fourier  coefficients  in  series  expansion 

for  N  (x,y);  see  equation  (B24) 
x 

Fourier  coefficients  in  series  expansions 

for  N  (x,0)  and  N  (x,b)  respectively;  see 
x  x 

equations  (B25)  and  (B26) 


_  j  ii 

V  Sm 


gm  C1  mit/a  >  em  C1 

thickness  when  plate  is  isotropic 


1,2,5,  or  4 

Fourier  coefficients  in  series  expansion 
;  see  equation  (B29) 

Fourier  coefficients  in  series  expansion 

for  N  ;  see  equation  (B27);  see  equation 
xy 

(B58)  for  value  of 

combinations  of  known  Fourier  coefficients, 
defined  by  equation  (B67) 

known  quantity  defined  by  equation  (B72) 


m,  n,  p,  q 


summation  indexes  (integers) 


upper  limit  on  m,  p,  and  q 
summation  index  (integer) 


upper  limit  on  n 
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N]_(y)>  N2(y)>  n?(x),  n^(x) 

external  running  tensions,  force  per 
unit  length;  see  figure  1 

N  ,  N  ,  N 
x7  y  *y 

plate  stress -resultants,  force  per 
unit  length;  see  figure  3 

P 

summation  index  (integer) 

Px(y),  P2(y),  Vx)’  Fh^ 

stiffener  cross-sectional  tensions; 
see  figure  3 

p1(o),p1(b),  p2(o),  P2(t>),\ 

P5(0),  P3(a),  P4(0),  P4(a )/ 

stiffener  end  loads;  see  figure  1 

P 

common  value  of  the  above  when  all  are 
equal;  used  in  numerical  example 

q 

summation  index  (integer) 

^(y)*  ^(y),  q5(x),  q^(x) 

external  shear- flow  loadings;  see  figure  1 

^  ^7  7  <C' 

Fourier  coefficients  in  series  expansions 
f°r  q1,  q_2,  q.y  q^  respectively;  see 

equations  (4)  and  (7) 

K’  Bn’  K"  ’  C 

combinations  of  known  Fourier  coefficients , 
defined  by  equations  (B 66)  and  (B67) 

s»,  a",  s"*  ,  a"" 
n7  n7  m  7  m 

Fourier  coefficients  in  series  expansions 
for  the  stiffener  cross-sectional  tensions; 
see  equations (BI7) 

t»,  t",  tw  ,  t"" 
n  n  m  7  m 

Fourier  coefficients  in  series  expansions 
for  the  derivatives  of  the  stiffener  cross- 
sectional  tensions;  see  equations  (B20) 
and  (B42) 

T 

mn 

Fourier  coefficients  in  series  expansion 
for  82ey/8x2  +  d2ex/$jr2  ;  see  equation  (6), 
also  equations  (8),  (9),  and  (10) 

rpf  mtt  mfff  m,,lt 

n7  V  m  7  rm 

Fourier  coefficients  in  series  expansions 
for  thermal-strain  discontinuities  between 
stiffeners  and  plate  edges;  see  equations 
(5)  and  (7) 

U,  V 

x  and  y  components  of  displacements  in  plate  j 

14 

IJT  UM 
mMN*  mMN 

v 

^Mmnp 


W1  w 

MNmp*  MNmp 

x' 


yr  y*1  y  ftt  ymm 

Mn*  Mn*  mii*  rJm 

y 

r,  yM 

Y 1  Yu 
Mn’  Wm 

71  7MI 

nto*  nm 


m 


a 


known  quantities  defined  by  equations  (B7I) 

plate  displacement  component  in  y-direction 

known  quantity  defined  by  equations  (B76) 
and  (B77) 

known  quantities  defined  by  equations  (B75) 

Cartesian  coordinates;  see  figure  1 

dummy  variable  representing  x 

known  quantities  defined  by  equations  (B73) 

Cartesian  coordinate;  see  figure  1 

dummy  variables  representing  y 

known  quantities  defined  by  equations  (B74) 

known  quantities  defined  by  equations  (B73") 

coefficient  of  thermal  expansion  of  plate 
and  stiffeners  in  numerical  example 


^nO*  ^mO*  ^mp 
b 

mpq 

An 

ex(x,y),  e  (x,y),  7xy(x,y) 
^(y),  e2(y),  e^x),  e^(x) 
e 


v 

^(y),  o*2 (y),  ffj(x),  o^(x) 
X 


KroneckerTs  delta;  unity  when  both  subscripts 
are  equal,  zero  otherwise 

1  for  m,  p  both  even  and  q  odd; 

1  for  m,  p  both  odd  and  q  even; 
zero  otherwise 

known  quantity  defined  by  equation  (B69) 
plate  total  strains 
stiffener  total  strains 

temperature  rise  of  plate  center  relative 
to  the  edges,  used  in  numerical  example 

Poisson’s  ratio  when  plate  is  isotropic 

prescribed  stiffener  stresses,  positive 
for  tension 

4ah/(jr2A);  area-ratio  parameter  used  in 
presenting  the  results  of  the  numerical 
examples  - 
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APPENDIX  B 


ANALYSIS  FOR  THE  CASE  OF  CONSTANT-AREA  STIFFENERS 


Basic  equations,  -  With  u(x,y)  and  v(x,y)  denoting  the  x-  and 
y-components  of  infinitesimal  displacement,,  the  strain-displacement 
relations  for  the  plate  are 

.  _  ^  Sv  Su  ,  Sv  f  s 

6x  "  ^  ey  =  ¥  7xy  =  ^  ^  (B1) 


Equations  (Bl)  imply  the  following  compatibility  condition  on  the  strains 


527 

xy 


The  plate  equilibrium  equations,  namely 


Sn 


x 


^r 


cN  SN 

7c-£  +  =  0 

qy  dx 


imply  the  existence  of  a  stress  function  F(x,y)  such  that 


(B2) 


(B3) 


Nx  -  S2F/^y2  N  =  52F/Sx2  N  =  -  S2f/ ck cy  (b4) 

y  xy 

Eliminating  the  strains  in  equation  (B2)  by  use  of  equations  (2),  and 
then  the  stresses  by  use  of  equations  (b4)  leads  to  the  following  form 
of  the  compatibility  condition,  in  which  account  is  already  taken  of  the 
equilibrium  and  stress-strain  relations: 


aV 

i? 


(C>,  - 


2C3) 


jtl_ .  +  c 

Sbc2^2  1 


a2e  a2e 


X 


ax2 


V 


=  0 


(B5) 


Considering  now  infinitesimal  lengths  of  the  stiffeners  as  free  bodies, 
and  utilizing  the  third  of  equations  (b4)  to  express  N  at  the  plate  edges 
•  4-  xy 

m  terms  of  F,  one  obtains  the  following  equilibrium  equations  governing 

the  longitudinal  variations  of  the  stiffener  cross-sectional  tensions: 
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dP-j/dy  -  (SaF/Sx^y)x=0  -  q1(y)  =  0 
dP2/dy  +  (S2F/Sx^)x=a  +  q2(y)  =  0 
dP^/dx  -  (daF/dx<V)y=0  -  q^(x)  =  0 


(B6) 


dP^/dx  +  (d^/dx^y)  +  q^(x)  =  0 

Integral  attachment  "between  the  stiffeners  and  the  plate  edges  implies 
equality  of  their  longitudinal  strains  and  leads  to  the  following 
additional  set  of  conditions  on  P^,  Pg,  Py  and  P^: 

P1(y)/(A1E1)  +  e1(y)  =  ey(0,y)  =  (ey  +  C2Ny  -  C5Nx)x=Q 


P2(y)/(A2E2)  +  e2(y)  =  6y(a,y)  =  (ey  + 


P^(x)/(A^EX)  +  e,(x)  =  eY(x,0)  =  (ex  +  -  C5Ny)y=0 


'  ^  y  3 


Vx)/(A4E4)  +  “  €Y(x/b)  "  (ev  +  "  C^Nv^v=' 


x  lx  3  y  y=b 


in  which  equations  (2)  have  been  used  to  obtain  the  right-hand  terms. 

The  assumption  of  negligible  bending  stiffness  for  the  stiffeners  permits 
the  substitutions  (Nx)x=0  =  N  (y),  etc.  to  be  made  in  these  terms.  Express 
ing  the  remaining  stress  resultants  in  terms  of  F,  one  obtains 


p-i  (y)  x 

* — e~  +  ■  ev(°>y)]  •  c2  +  c3Ni^y^  =  0 

A1  E1  1  y  2  &2  x=0  y  1 


♦  t.*w  -  +  °*(y) ' 0 


P,  (x) 

-p—T  +  [e,(x)  -  e  (x,0)]  -  C  (— )  +  C  N  (x)  =  0 

A3  E3  5  Sy2  y=°  0  ? 


T-  E~  +  [eJx)  -  e  (x,b)]  -  C  +  C  N^(x)  =  0 

\  E4  4  x  1  5y2  y=b  5  4 
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The  problem  can  now  be  stated  essentially  as  follows:  to  solve 
equations  (B5),  (B 6),  and  (BT)  for  F,  F±,  Pg,  P  and  P^,  subject  to 

boundary  conditions  arising  from  the  prescribed  forces  at  the  stiffener 
ends  and  the  prescribed  distributed  loadings  Ng,  N^,  and  N^.  In 

the  following  sections  a  formal  solution  to  this  problem  will  be  obtained 
in  terms  of  Fourier  series. 


Boundary  values  of  F.  -  The  fact  that  the  distributed  loadings  N  , 

W.and  W4  are  transmitted  directly  to  the  plate  means  that  the  second 
derivative  of  F  in  the  direction  of  the  edge  is  known.  Therefore  two  inte¬ 
grations  will  give  the  variation  of  F  along  each  edge  in  terms  of  the  unknown 
corner  values  and  the  known  N^,  Ng,  N^,  .  For  example  (using  subscripts  on 

F  now  for  convenience  to  denote  partial  differentiation), 


FyyCo.y)  =  Nx(y) 


Therefore 


Fy(0,y)  =  Fy(0,0)  +  J  N^(y'  )dy* 


Therefore 


y  y' 


F(°,y)  =  F(0,0)  +  yFy(0,0)  +JJ  \  (y"  )dy"dy 


Substitution  of  y  =  b  in  equation  (b8)  gives 


b  y' 


Fy(°,°)  =  £  [  F(0,b)  -  F(0,0)  -  J  J  N1(y")dy"dy'  ] 


which  result,  substituted  back  into  equation  (B8),  gives 

by'  yy' 

F(0,y)  =  F(0,0)  +£  [F(0,b)  -  F(0,0)  -  J  J  N1(y")dy,,dy'  ]+J 'J  N1(y")dy"dy 


Thus  the  variation  of  F  along  the  edge  x  =  0  has  been  determined  to  within 
two  constants,  F(0,0)  and  F(0,b).  Replacing  N^(y”)  by  its  series  expansion, 

equation  (3),  and  carrying  out  the  integrations  indicated  in  equation  (B9) 
give 


F(0,y)  =  F(0,0)  +  £  [F(0,b)  -  F(0,0)]  -  L  B'  (— ^sin  ^ 


n=l  n  nn 


(BIO) 
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going  through  a  similar  procedure  for  each  of  the  remaining  edges,  one 
obtains 

F(a,y)  =  F(a,0)  +  £  [F(a,b)  -  F(a,0)]  -  ng1  (^)  sin  ^  (Bll) 

F(x,0)  =  F(0,0)  +  f  [F(a,0)  -  F(0,0)]-  B”'  (^)  sin  ^  (B12) 

M  2 

F(x,h)  =  F(0,h)  +  |  [F(a,h)  -  F(0,b)]  -  ^  B^"  (^)  sin  2HHS  (Blj) 


For  later  use  it  will  he  necessary  to  have  Fourier  expansions  of  the 


'  in  the 

form 

F(0,y) 

ii 

a  * 
n 

sin 

(njty/b ) 

F(a,y) 

N 

=  n=l 

an 

n 

sin 

(nity/b) 

F(x,0) 

M 

= 

m~l 

a,!l 

m 

r  sin 

(mnx/a) 

F(x,b) 

1! 

0' 

1IMS 

a,M 

m 

'  sin 

(mnx/a) 

(bi4) 


Evaluating  the  coefficients  in  these  series  through  the  formulas 
rb 

a'  =  (2/b)  /  F(0,y)  sin  (nny/b)dy,  etc.,  with  F(0,y),  etc.  replaced 

0 

by  the  right-hand  sides  of  equations  (BIO)  to  (B13),  one  obtains 

an  ■  I;  (F(0>0)  -  (-1)n  F(0’b)I  • 

^  =  h  [F(a’0)  -  (-1)n  F(a'b)!  -  B» 

a”'  =  h  [F(°>0)  “  (_l)in  " 

a""  =  —  [F(0,b)  -  (-l)m  F(a,b)]  - 
m  .  mit 


(— ) 

ttitt 


B”  1 

m 


Bf 


m 


(B15) 
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Series  assumptions  for  F(x,y)  and  P^  Pg,  P^,  P^.  -  In  the  interior 

the  plate  (i.e.  in  the  region  0  <  x  <  a,  0  <  y  <  b,  excluding  the  edges 
x  _  a  and  y  =  0,  b),  the  stress  function  F(x,y)  will  be  assumed  to  be 
representable  by  the  double  Fourier  series 

.  M  N  . 

F(x>y)  =  n|1  a^  sin  (mjcx/a)  sin  (nny/b)  (Bl6) 


with  as  yet  unknown  coefficients.  Equation  (Bl6)  is,  of  course,  not  -valid 
at  the  edges;  however  there  the  values  of  F  are  already  expressed  in  series 
form  by  equations  (Bl4)  and  (B15).  Similarly  the  stiffener  forces  will  be 
assumed  in  the  form 


Pq(y)  = 

N 

:  nil 

s 1 

n 

sin 

(noty/b) 

P2(y)  = 

N 

:  nil 

s" 

n 

sin 

(n:rty/b ) 

P3(x)  = 

M 

m=l 

s,n 

m 

f  sin 

(mirx/a) 

P400  = 

M 

mil 

s,m 

m 

sin 

(nutx/a) 

(B17) 


for  all  cross  sections  except  the  end  cross  sections.  At  the  end  cross 
sections  the  stiffener  forces  are  already  known  from  the  given  loading 
(see  fig.  1). 


The  coefficients  in  the  series  in  equations  (Bl6)  and  (BI7)  are 
lated  to  the  left-hand  sides  through  the  usual  formulas 


a  b 


mn 


=  —  r 

ab  J 


F(x,y)  sin  222.  sin  22£  dydx 
8.  D 


s  *  =  § 


Pq(y)  sin  25£  dy  f  etc. 


(B18) 

(B19) 


Series  for  the  derivatives  of  F(x,y)  and  P^,  P£,  P^,  P^.  -  The 

derivatives  appearing  in  equations  (bV)  to  .(B7)  will  also  be"  assumed 
expressible  in  series  as  follows 
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d.P1/dy  =  n2^  V  cos  (mty/b)  (0  <  y  <  b) 


dPg/dy  =  n|0  t"  cos  (rnty/b)  (0  <  y  <  b) 


jyL 

dP,/dx  =  2L.  t,,f  cos  (irurx/a)  (0  <  x  <  a) 

y  m=0  m  '  -  — 


1V1 

dP,  /dx  =  Z.  tM"  cos  (mnx/a)  (0  <  x  <  a) 
4 7  m=0  m  '  —  - 


(B20) 


M  N 

(.  =  8%/^  =  n£l  ngl  V  sin 


(mitx/a)  sin  (nity/b) 


0  <.  x  <  a 


0  <  y  <  b 


(B21) 


(^y)x_0=  ( B2F/ck2) 


_  =  Z.  c'  sin  (nity/b) 
x=0  n=l  n 


(0  <  y  <  b)  (B22) 


(N  )  =  (82F/ck2)  =  Z  c"  sin  (nrty/b) 

v  y  x=a  v  '  x=a  n=l  n 


(0  <  y  <  b)  (B25) 


M  N  <  x  <  a1 

N  =  d2F/cy2  =  Sm„  sin  (nutx/a)  sin  (nity/b) 

x  J  m=l  n=l  inn  lo  <  y  <  b. 


(B24) 


(Nx)y=0=  ^2F/^2)y=0  =  mil  g4  Sln  (lmtx/a) 


(0  <  x  <  a)  (B25) 


(Nx}y=b=  (a2F/^2)y=15  =  mil  gm  sin  (mrtx/a) 


(0  <  x  <  a)  (B26) 


-  »/**  ■  j0  x  j"»  °oB- (mxM  °°s  (n’y/b)  !  f  4 


(B27) 


d^F/dx*  =  m51  sin  (nutx/a)  sin  (nity/b) 


M  N 


0  <  x  <  a 


0  <  y  <  b 


(B28) 


c^F/dy*  =  Z,  Z,  i  sin  (mjtx/a)  sin  (nity/b) 

»  ^  in  — 1  -n— 1  mn 


M  N 


m=l  n=l  mn 


0  <  x  <  a 


,0  <  y  <  bj 


(B29) 


M  N 


5  F/dx2^2  =  E_  L  V  sin  (imtx/a)  sin  (nity/b) 
7  m=l  n=l  mn 


0  <  x  <  a 
0  <  y  <  b 


(B50) 
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where 


22 


in  the  series  for  the  basic  quantities  (equations  (Bl6)  and  (B17)).  The 
former  can  be  expressed  in  terms  of  the  latter  by  means  of  integrations 
by  parts  in  the  right-hand  sides  of  equations  (B3l)  to  (B4l)*.  For 
example,  from  equations  (B3l) 


2-6 


[P^(b)  cos  mt 


P,  (0)  +  /  p  (y)  sin  s 2L  dy ] 


2-6 


P,(b)  cos  ntt  -  Pn(0)]  +  y 


2-6 


[P0(b)  cos  mt  -  Po(0)]  +  s'^ 


(B42) 


2-6 


[P5(a)  cos  mn  -  P^O)]  +  f-  s"* 


2-5 


m  a 


[  P),  (a)  cos  irnt  -  P.(O)]  +  7-  s"" 


Similarly  two  partial  integrations  with  respect  to  x  in  equation  (B32) 


c  =  —  -  [a'  -  (-1)“  a"]  -  (2L)  a 
ran  a  a  L  n  v  7  nJ  va  ran 


m  _  /-rait 


(b43) 


Two  with  respect  to  y  in  equation  (B35)  give 


=  ^  |  [a”'  -  (-l)n  a""]  -  (SS) 


’ran  b  b  m 


(B 44) 


In  equation  (B38)  partial  integration  with  respect  to  x,  followed  by 
partial  integration  with  respect  to  y  in  both  of  the  resulting  terms,  gives 

J  =  -  &m°^2  &no)-[  (-i)m+n  F(a,b)-(-l)m  F(a,0)-(-l)n  P(0,b)  +  F(0,0)] 
ran  ab 


nut  f^nO  f  ( _x)n  ,,,,  _  a,M]  +  mt  f^mO  f  „  _  a 

a  b  L  v  J  in  m  b  a  n 


"  -  a']  +  -^a  (B45) 

n  n  a  b  ran 


Such  a  technique  was  employed  for  plate  bending  problems  by  A.  E.  Green 
(reference  l),  who  ascribes  its  earlier  use  to  S.  Goldstein. 


in  which  single -valuedness  of  F  at  the  corners  has  been  assumed.  Proceed¬ 
ing  in  a  similar  fashion  with  the  right-hand  sides  of  equations  (B39)  to 
(B4l),  one  obtains 


■Z  L 

e  =  2HL  ^  r  c  1  _  /  n  iti  /nut\  2  r  .  /  n  vui  ,,,  /mjr\ 

emn  a  a  Lcn  '  V  "  a  [an  “  (_l)  an]  +  V 


(b46) 


5 

i  =  ££  —  [g1  -  (~l)n  g"l  -  C— )  —  ra'"  -  (-1  'i11  n""l  + 

mn  1)  t  lsm  K  XJ  SmJ  ^  '  b  lam  ^  1}  a™  J  + 


m  J  'b  '  mn 


(B4T) 


pim  "  ST  T  f  1  *(«>!>)  -  (-1)”  F(a,0)  -  (-1)“  F(0,b)  +  F(0,0)] 


+  I  (?)  ft  (-1)”  C-»r)+  1(f)'  =  [  C-Dm  <-a;]+(?)a(S?) 


2 - 2 


m  m  a'b  '  a 


n  n  'a  'b  '  “mn 


(B48) 


Substitution  of  series  into  the  basic  equations.  -  Through  equations 
(B42)  to  (B48)  all  the  unknown  coefficients  in  the  derivative  series  are 

expressed  in  terms  of  the  basic  unknowns  a  :  cf ,  c" ,  s' .  g":  s' .  s"  s,u 
»  %  1  .  .  mn  n*  n}  Bm>  V  n'  V  m  ' 

sm  *  F(a^O),  F(0,b),  and  F(0,0).  Relationships  among  these  unknowns 

will  now  be  obtained  by  substituting  the  assumed  series  into  the  basic  equa¬ 
tions  (B5),  (B6),  and  (B7) . 

Considering  first  equation  (B5),  substituting  into  it  the  series 

expansions  from  equations  (6)  and  (B28)  to  (B30),  and  eliminating  e  ,  i  , 

mn'  mn' 

and  p  through  equations  (B4b)  to  (B48),  one  obtains 


n  Jmjt  2  r  , 


K  -  (-1)”  -  (f  )3  I  K  -  (-1)”  O  *  (f  ^ 


+  (°4  -  *-.)  {^7X  [(-l)m+nF(a,b)-(-l)"V(a,0)-(-l)n  F(0,b)  +  F(0,0)] 


+  f  Off  f  [  (-D”  a""  -  a'"J  +  |  &.f  ?I  (-1)"  a"-a']+(=)'(S) 


2  „  2 


m  m  a  xb  7  a 


n  nJ  'a  '  xb  '  inn 


*  ci  f  K  ■  (-1)”  -  <f )?  I K'  -  (-«" »;")  *  O  U  p  V  -  o 


Solving  this  equation  for  a^  and  eliminating  at,  a^,  ,  and  a""  through 

equations  (B15),  one  obtains 


2k 


mn  2 

mn  jt 


—  [ (-l)m+n  F(a,b)  -  (-l)m  F(a,0)  -  (-l)n  F(0,b)  +  F(0,0)] 


_1_  L  +  £  mil  [c.  _  (_i)m  c"]  C.  +  |  ^  [g*  -  (-l)n  g"l  C 
E  1  mn  a  a  L  n  v  '  nJ  2  b  b  °m  m  1 

mn  L 


I  T  ^)P[Bn  ~  ^ C2  +  O  -  2C5)] 


- 1  r  &>  ^  ci +  (f}  (c*  -  2C3)] 


(b49) 


where 


/mt\2  .  „  /nir 


V  ■  c2  O  *  <C4  -  “,)  <?>  <?  >  +  C1  ‘w' 


(B50) 


Thus  through  the  compatibility  equation  the  unknown  a^  have  been  expressed 
in  terms  of  a  smaller  class  of  unknowns ^  namely  the  c^,  c”,  and  g^. 

Turning  now  to  the  stiffener  equilibrium  equations  (b6),  substituting 
the  series  from  equations  (B20),  (B27)j  and  (4),  and  utilizing  equations 
(B42),  one  obtains  the  relationships 


2-6 


“°  [(-!)”  PjCb)  -  PitO))  +  V  Bn  -  %  '  Jo  >m  ~  0  fr-0-1* 


(B51) 


2-8 


n0  [(-l)n  P2(b)  -  P2(0)]  +TSn  +  ^  +  mSo  3mn(-l)ia=  ° 

(B52) 


2-6  N 

— ^  [(-i)"  p,(«)  -  f5(°)i  f »;n-  <c-  .fc  ■ 0 


(B53) 


„5-°  [(-l)m  P,,(a)  -  P,,  (0)3  +  ^  s"”+  Q""  +  jmn(-l)n  =  0  (m=0,l,.-M) 


(B54) 


The  n  =  0  and  m  =  0  equations  of  this  group  will  be  studied  first.  From 
equation  (b45),  in  conjunction  with  (B15)  and  (b49)j  it  is  noted  that 
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^00  ab  "  ^(a,0)  -  F(0,b)  +  F(a,b)] 


(B55) 


'On  a  m  (BA  "  Bn)  for  n  ^  0 


(B5  6) 


mO  b  mjT  x  m 


t  -  C>  for  m  ^  0 


Therefore  equations  (B51)  and  (B52)  for  n  •  0,  (B53)  and  (B5W  for 
m  =  0  give 


Px(b)  -  P1(0)  -  bQ^  -  bj00  +  J^a/m*)^"  -  B™  )  =  0 


(B57a) 


M 

f2(b)  -  P2[0>  +  bQj  +  bj00  -  ^(a/mnHB-  -  B;")(-l)m  = 


(B57b) 


P3(a)  "  P3(0)  -  a%"  -  aj00  +  nt  (b/n*)(Bn  "  ^ )  = 


(B57c) 


I\T 

Va)  -  P4(0)  +  a<^"  +  ajQ0  -  nZ1  (b/n«)(B"  -  B> )  (-l)n  =  0 


(B57d) 


Three  of  equations  (B57)  are  redundant  if  the  applied  loading  constitutes 
an  equilibrium  system  for  the  structure  as  a  whole:  It  is  easily  verified 
that  the  equation  of  equilibrium  of  forces  in  the  x-direction  is  the  same 
as  that  obtained  by  adding  equations  (c)  and  (d),  the  equation  of  equilibrium 
or  the  y-direction  the  same  as  the  sum  of  equations  (a) .and  (b),  and  the 

fsSrPH'  M?ntS  ^OUt  the  center  of  the  plate  equivalent  to 

(a/2)[eq.(a)  -  eq.(b)]  .+  (b/2)[eq. (d)-eq.  (c)] ,  Hence  any  three  of  equations 
^B57;  may  be  eliminated  from  further  consideration.  The  one  remaining  eaua- 
lon  serves  to  establish  the  value  of  j...  Selecting  equation  (a)  as  this 
remaining  one  gives  00 


(B58) 


Three  alternate  expressions  for  j  can  be  obtained  from  equations  (b), 
^c;,  and  {&).'  For  the  sake  of  symmetry  one  can  sum  all  four  expressions 
.r.J00  and"  divid-e  ,by  four  ln  order  to  obtain  still  another  formula  for 
Jqq)  while  esthetically  more  satisfying,  this  formula  is  more  complicated 
than  equation  (B58) . 

Finally,  equations  (B7),  with  the  various  terms  replaced  by  their 
series  expansions  from  equations  (3),  (5),  (B22),  (B23),  (B25),  and  (B26), 
give 
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°y<Vi> +  t;  -  c2=; +  c3b;  ■ 


s"/(a2e2)  +  t;  -  c2c;  +  c  b;  -  0 


s»./(A5e3)  +  T-  -  CA  +  CjB”'  -  0 


(B59) 


C/(aA)  +  t;-  -  cA  *  o3b;"  -  o 

In  obtaining  equations  (B59),  the  assumption  of  uniform  stiffeners  was  used 
for  the  first  time. 


Reduction  in  the  number  of  simultaneous  equations.  -  Equations  (B5l) 
to  (B54),  with  the  n  =  0  and  m  =  0  equations  excluded,  can  he  used  to  obtain 
a  system  of  simultaneous  equations  in  which  the  c^,  c^,  and  are  the 

only  unknowns.  This  is  accomplished  by  eliminating  s^,  s^,  s^f  y  an(^-  sm 

the  aid  of  equations  (B59) >  an(i  means  of  the  following  expression 

=  .  f  +  2  mr  [c,  _  (_±)m  ^  c  +  £  M  [g--  .  (-l)ng^]C  } 

Jmn  ab  E  |  mn  a  a  n  "  ° k  *h  m  ™ 


n  2  b  b  m 


+  —  -p  (— )  [b"»  -  (-l)n  B""]  C0  + 
E  |b  '  a  '  L  m  v  '  mJ2 
mn  L 


f  (f)5[B’-(-l)V]Ci}  (b60> 


which  is  obtained  from  equations  (b45),  (B15),  and  (B49) •  The  resulting 
system  of  simultaneous  equations  is 

c’  [A  E  +  —  !  Waif]  .  c"  [2  XL  W»)‘(-1)!] 

A  1  A  »  -1  E™;  A  »  "-1  W. 


O  2  M  8 ’  -  (-1)  g^  .  . 

=  R*  -  |  (^)  E,  ~  g - —  (n=l,2,...,N) 

n  h  vb  m  1  V 


(b6i) 


r2  M  +  ;  E  +  |  «  £hM 

J  n  L  2  2  a  m=l  E 


c«  r-  z  >  7 ■  ■ 
Cn  a  m=l  E 


K  * 1  I  (-1)”  ~  2  - 


(B62) 
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'  [A  E  +  I  I,  -(nrtA)2i  r2  J  (n« /b)a  (-l)n, 

m  5  3  b  n=L  E  J  "  sm  '■‘b  n=l  E  J 


=  Rm  _  2  /mrtv2  |  °n  "  Cn  ,  * 

m  a  '  n=l  E  (m-1,2, . .  .,M) 

mn 


(B63) 


’  [f  2,  (nit/b)g  (-1) 


lr  U  v  V--W  i  ,  ta  w  x  c  v  ■ 

m  B  n=1  E  J  gm  tA4E4  +  b  n=l  E 


(nit/'b) 


=  r"«  +  £  (— )2  Z  (-I)n  °n  "  ^~1') 

m  a  v  a  '  n=l  v  '  E 

mn 


(m=l,2, . 


(B64) 


where 


a  CnC2  n;t/b  '  cn  "  CnC2  nit/b  '  =  gmCl  mit/a  '  Sm  =  gmCl  mir/a  (b65) 


and  Rn’  RnJ  Rm"-’  Rm"  are  the  follovinS  combinations  of  known  thermal  and 
loading  quantities: 


x  2 


a  =  ^  +  p.Cb)]  +  ±  (b;-b;)  +  aA  sg(c  BN*y  -  J, 


Rn  =X  +F[P2^-(-1)n  p2(b)l  - 


(B’-B")  +  APE?  t~(C,B"+Tm)  +  Zn  (-1)“k 
anjt  n  n  d  d  h  '  5  n  n'  m=l  mn 

(b66) 


E»'  -  V  *  !'e3(0)-(-i )"  P3(«)J  *  £  (B"'-B"")+  Vj  =(c3b;..  +C)- 


. t* tt  _  „l 


in.  2 


-s»  +  ^  V°M-1)  Va»  -  H(c3b;"+t;")  * 


'-■r{vVE  0\  [B"'-(-Dn  b;-i-  |(f!)sc1[B^-(-i)V)|  (b67) 


Equations  (B6l)  to  (B65)  can  be  solved  simultaneously  for  the  c ' ,  c" , 


g^>  and  g".  With  these  known,  equations  (B59)  will  furnish  the 


values  of 


Sn’  Sn>  Sm  ’’  Sm" >  and  e<luati°ns  (b6o),  (B58),  and  (B56)  the  values  of  the 


28 


j  .  Equations  (B17)  will  then  give  the  stiffener  stresses,  equations  (B2l) 
to  (B27)  the  plate  stresses . 

Further  reduction  in  the  number  of  simultaneous  equations.  -  Equations 

(B6l)  and  (B62),  written  for  the  same  value  of  n,  can  be  solved  for  each 

and  c"  in  terms  of  all  the  g’  and  g",  with  the  result 
n  mm 

2  s'  -  (-l)n  g" 

1  J,.T  ,2  M  (nut/a)  1r  ,  2  rmtx2  M  fm _  j  .V 

cn  "  A  |[A2E2  a  m=l  E  ^Rn  '  h  ^b  '  m=l  E 

n  ^  mn  LmL 


f2  M  (mjt /a)  (-!?”] [R»  +  i  (M)2  £  (_i)m  °m 

L  a  m^L  E  n  b  b  m=l 

mn 


K  -  <-1)n  K 


E 


mn 


?■  .  i  {[AA  ♦  f  l  *  I  1 

n  An  [l  1  1  a  m=l  E^  n  b  h  m-i. 


(B68) 


+  .2  M  (W.]i(-ir,r,  _  2  (SL)*  M  -  (-1)n  ? 

+  ^a  m=l  E  n  b  ^b  ‘  m=l  E 

1 


m- 


mn 


mn 


where 


r  A  „  2  ^  (nrn/a)  1Pa  n  .2  ^  (mrt/a)  ■.  r2  ^  (mrt/a)  (-1)  ^ 

^  -  "Vi  W;Jl  E]m  1 

Using  equations  (b68)  to  eliminate  the  c*  and  c"  in  equations  (B63)  and 
(b64),  one  obtains  the  following  system  of  simultaneous  equations 
involving  only  the  g1  and  g"  as  unknowns: 


?  obtains  the  following  system  c 
only  the  g^  and  g^  as  unknowns: 

Si  <wwmp  -  smp  *&>  +  i;  (»*mp  +  5„p  ]  -  - 

M  r  _  _ 

pil  [  gp  ^WMKmp  +  5mp  YNm^  +  gp  ^WMNmp 


M 

pil 


U 


mMN 


(B70) 


8mp  XNm^ ] 


-  u" 

mMN 


(m  =  1,  2,  M) 
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where 


IT'  =  rim  _  2 

mMN  m  a  'a  '  n=l  E  A 


mn  n 


U".„T  =  R""  +  -  (— )2  l  (_l)n  hin 


mMN  m  a 


n=l  E  A 
mn  n 


(B71) 


htn  *Mn  En  +  YMn  Rn  "  [X4i  Rn  +  Y4r  Rn  ] 


(B72) 


^Mn  AlEi  +  a  p?i 


2  M  (^/n\z 


p=l  E 


pn 


=  A„E„  +  £  l£2Zalf 


2  2  a  p=l  E 


pn 


x;  m  =  a  e  +|  I  InTr/b) 

Mm  3  5b  n=l  E 


mn 


(B73) 


y""  -  &  t?  j_  2  N  (rut/b) 

^m  -  W  +  b  nk.  E-1 

mn 


Y'  =  —  2  (P"A)aM)p 

Mn  a  p=l  E 

pn 

Y"  =  |  I  _(mt/b)a(-l)n 
Nm  b  n=l  E 

mn 


(B74) 


W'  =  iL  ^  (mt/b)2 

i®imp  ab  '‘a  '  n=L  E  E  A  V1 


mn  pn  An  ^ 


W"  =  _  iL  (EE)2  v  (-l)n(nn/b)a 
MNmp  ab  '  a  '  n=l  E^  E^  An  VMmnp 


(B75) 
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w  A2E2  +  AiV-1>”+p  + ;  Jl  v 


(B7  6) 


[1  -  (-l)p+<1  -  +  (-l)m+p] 


1  for  m  and  p  both 

even,  q  odd  /B77\ 

1  for  m  and  p  both  [  ' 

odd,  q  even 
0  otherwise 


The  advantage  of  this  preceding  step  is  evident:  Whereas  the  original 
system,  equations  (B6l)  to  (B 64),  requires  the  solution  of  2N  +  2M  simul¬ 
taneous  equations,  the  reduced  system,  equations  (B70),  contains  only  2M 
simultaneous  equations.  Thus  N  may  be  taken  arbitrarily  large  without 
increasing  the  number  of  simultaneous  equations  that  have  to  be  solved. 

Special  case:  symmetry  about  y  -  b/2.  -  If  the  structure,  loading, 
and  thermal  strains  are  symmetrical  about  the  line  y  =  b/2,  then  a  corre¬ 
sponding  symmetry  obtains  in  the  stress  function  F  and  in  the  plate  and 
stiffener  stresses.  Consequently  one  may  set 


c !  =  c  —  0 
n  n 


=  si 


for  n  even 


for  all  m 


(278) 


In  place  of  equations  (B6l)  to  (b64)  the  following  system  results: 


c'  [A  E  +  -  Z.  Wa)2]  _  c"  [£  §  (mrt/a)a(-.ll..; 

n  L  1  1  am=lE  •  nLa  m=l  E _ 

mn  mn 


K  -  E  <T>“  Jl 


(b6i*  ) 


(n  =1,  3,  N) 


o'  [S  1  +  e;  [A2E2  +  f  X¥^ 

n  a  m=!  n  2  2  a  m-1 


=  R*  +  £■  (^)2  1  (-i)m  (eVE  ) 

n  b  b  '  m=l  m  11111 


(b62‘  ) 


(n  —  1, 3>  •  •  •  >  E) 
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n,x 

mn 

=  R"' 

_  S  /m«\2  ^ 

a  V  a '  n=l,3, . . 

c'  -  (-l)m  c" 
n  n 

m 

E 

mn 

(m 

=  2^  •  • « 9  m) 

(b Gy) 


where  N  is  now  restricted  to  odd  integers. 

obtainedrSP°ndinS;i;7'  “  PlaCS  °f  equations  (B?°)  the  following  system  is 

M  _ 

P=1  Sp  (2WMNmp  +  5mp  YNm  '  8mp  ®  ^  (B70') 

(m  =  1,2,...,  M) 

Symb°lS  are  de fined  as  before  through  equations  (B7l) 

with  resi  '  bUt  Wlth  thS  understanding  that  wherever  a  summation  ' 
j|Jh_  ®  pect  n  appears  ln  these  definitions  n  shall  now  take  on  onlv 
d  values.  The  number  of  equations  in  the  system  (B70')  is  M  regard¬ 
less  of  the  value  selected  for  N.  ’  g  d 

11-1  LL-fl”1!.1  Case;  Bynmetry  about  y  »  b/2  and  x  =  a/P.  .  A  further 
loading  aJd  J5®  aUfb6+  °!  slmultaneous  equation  results  if  the  structure, 

y .  atout  both  “nterii'»8’ 


c  r 
n 

=  c" 
n 

=  0 

for 

n 

even 

g» 

=  0 

for 

m 

even 

I  0 

=  c" 
n 

for 

n 

odd 

gf 

=  I" 

°m 

for 

m 

odd 

(B79) 


Equations  (B6l>)  to  (B631)  are  replaced  by 
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k  M 
c 1  [  AnE  +  - 
nil  a  m=l,3> • •  E 


(irm/a) 


mn 


k  - 1  <f>UL.  <</E»> 

(n=i>3>  •  •  • )  n) 


(b6i" ) 


rA  E  +  -  I  (mM 
®m  3  3  13  n=l,3>  •  •  E 


]  =  R" ’  -  i  (E)2  I  (c'/E  )  (B63") 

J  m  a  v  a  '  n=173,..  n'  mn7 


mn 


(m-1,  3j  •  *  • )  M) 


where  M  and  N  are  now  restricted  to  odd  integers. 

Eliminating  the  cP  in  equation  (B63")  "by  means  of  equations  (B6lrt) 

gives  the  following  system  of  simultaneous  equations  in  the  g^  alone, 
which  takes  the  place  of  equations  (B7O  V) : 


M 

p=l^3, • • 


11  #  <=> 


p  •  ah  a 


tV’ib1'1  *  8mp 

-  a  <f  >2  n-lk-  Rn  "A'l 

(m=l,3»  •  •  •>  M) 


7  It  I 

Tim 


} 


where 


k  M 

=  A1E1  +  a  m=l?3, 
4  N 

Z'"  =  A_E_  +7-  ,2, 

Nm  3  3  n=l,3* 


(mjt/a)g 

E 

mn 


0W^)g 

E 

mn 


(B70") 


(B73") 


Regardless  of  how  large  a  value  is  assigned  to  N>  there  are  only  (M+l)/2 
simultaneous  equations  in  the  system  (B70"). 


33 


cass  of  large  stiffener  areas#  -  The  case  in  which  the 
stiffener  cross-sectional  areas  are  large  compared  to  the  plate  cross' 
sectional  area  is  of  practical  and  theoretical  interest*  In  order  to 
study  this  case,  let  it  be  assumed  that  AgEg,  and  A^  all 

approach  infinity  while  maintaining  constant  ratios  with  each  other. 
Then  equations  (B6l)  to  (b64)  degenerate  to 


°n  -  b;/<aa> 


■  r;"/(a3e3> 


g"  -  (R""/ (A,  Ei  ) 


(b8o) 


Thus  for  this  limiting  case  it  is  unnecessary  to  solve  simultaneous 
equations* 


Equations  (b8o)  represent  a  first-order  perturbation  solution  of 
equations  (B61)  to  (B64),  in  which  the  coefficients  of  R^/(A  Ep, 

RnAA2E2^  etc‘  are  expanded  in  series  of  powers  of  (a3  )_1  and 

only  the  zeroth  power  retained.  If  terms  of  the  zeroth  and  first  degree 
are  retained,  the  following  more  accurate  results  are  arrived  at: 


R' 


2  M  A 
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E  AnE, 
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2  /M)2  “  - — -  fJH _ 'i  x  2  ,n*w  ,\n  “  1  A’  s 

b  U  ^  m?l  AA  ^  b  (b  M  1)  m§!  Emn  AA  ^ 
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+  £  (— )2  i,  — (JL_)  _  2  ffi2,  ,n  M  (-l)m  A"  % 

b  b  rA  Em  A2E2  A^E^  b  (b}  (  }  mil  A^  A^E^ 


2/  n  \n 


«.  =  _J1_  1.2  N  Wh)a  +  Rrn"  2  N  (nTt/b)5^-!)1 
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(B84) 


It  is  seen  that  a  solution  to  this  degree  of  approximation  still  retains 
the  simplicity  of  equations  (b80) ,  in  the  sense  that  it  is  unnecessary  to 
solve  simultaneous  equations  to  determine  the  g^  and  g”*  The  coef¬ 

ficients  of  R'/(AnE  ),  etc.  in  equations  (B8l)  to  (B84)  are  correct  to  terms 
n  -L  i 

of  the  first  degree  in  (a3E  -A  E  )“  .  To  obtain  expressions  of  greater 
accuracy  it  would  once  again^be  necessary  to  solve  simultaneous  equations. 
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APPENDIX  C 


ANALYSIS  PX)R  THE  CASE  OP1  CONSTANT-STRESS 
(OR  PRESCRIBED-STRESS)  STIFFENERS 


h  stiffener  cross  sections  are  varied  so  as  to  produce  in  each 
tiffener  a  constant  (or  non-constant  but  prescribed)  stress  everywhere 
along  its  length,  then  the  determination  of  the  plate  stresses  becomes 
?n  ^mple>  for  ln  that  case  the  normal  stress  in  the  adjacent  plate 

quentS'the  c-  ^  ^  Stlffener essentially  given, 'and  conse- 

known  fThis11^  gm'.Srr,  ar\  e^atlons  (B22h  (B23),  (B25),  and  (B26)  are 
is  is  in  contrast  to  the  constant-area  case,  in  which  the  c'  etc 
were  unknown.)  With  these  quantities  known,  equations  (b6o),  (B58),  n(B56)* 

Sn  en  ^  ^  (B24),  and  (B27)  will  furnish  the  remaining  plate  stresses. 

To  carry  out  the  foregoing  analysis  in  detail,  let  o^y),  ^(y),  a  (x), 

(^(x)  denote  the  prescribed,  perhaps  constant,  values  of  the  stresses  (P./A.) 

in  the  stiffeners  at  x  =  0,  x  =  a,  y  =  0,  and  y  =  b  respectively.  Then  the1 
trains  m  these  stiffeners  are  known  from  equations  (l).  Equating  these  to 
the  corresponding  plate  strains  adjacent  to  the  stiffener,  equations  (2),  gives 


+  (ct1/E1)  =  e  (0,y)  +  C  N  (0,y)  -  C  N  (y) 


2  +  (o‘2/E2)  =  ey(a>y)  +  C2Ny(a,y)  -  C  N2(y) 


+  “  ex(x*°)  +  CpV^o)  -  C  N  (x) 

e4  +  (W  “  ex(x^b)  +  °A(X^)  -  C,N,(x) 


whence 


N  (0,y)  =  (2-1) 


dxs  X=0  ; 


c;  fei(y)  -  e  (0,y)  +  — -  +  C  N  (y)] 


(B3?w££!  tM0Wn  f  St“\iand  SMe"  °f  e'p:iations  (C2)  into  equations 
w  /  ,,  ’  '  ^  ^ ’  an<b  an<l  carrying  out  the  integrations,  one 

ms  the  values  of  the  c£,  c£,  g^,  and  g£,  from  which  the  plate  stresses 

can  be  found  as  indicated  in  the  previous  paragraph . 
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In  order  to  determine  the  variations  of  stiffener  cross-sectional  area 
required  to  produce  the  prescribed  stiffener  stresses,  the  relationships 


Ai(y)  =  P1(y)/cr1(y) 


A2(y)  =  P2(y)/<J2(y) 


Aj(x)  =  Pj(x)/cr5(x)' 


A^(x)  =  Pj+(x)/ci)+(x) 


may  be  used,  with  P^(y)  etc.  evaluated  by  integrating  the  running  loads 
on  the  stiffeners,  starting  at  one  end,  as  shown  below. 


Px(y) 


J  «/ 

Px(0)  +  J  g1(y)  dy'  -  J  Nxy(0,y)  dy 


(o)  +  J  q1(y')  dy-  +  j  (g^)  dy 


p  n  M  N 

=  P1(0)  +  J  q^y*)  dy-  +  J  mZ0  nZ0  jmn  COS  dy 


p  M  N  -u 

Pl(0)  +  f  %(y,)  dy’  +  mi0  (yJmO  +  n=l  m  sin  T) 


Similarly 


P2(y)  =  P2(0)  -  J  q2(y’)  dy  *  -  (-l)m(yjm0  +  n|_  ^  V  sin  ^L)  (C5) 


P  (x)  =  P  (0)  +  f  q (x!)  dx*  +  Z  (xj  +  Z  —  j  sin  22£) 
3V  y  3W  /  y  n=0  v  °0n  m=l  mrt  °mn  a  J 


\  (x)  =  Pj  (o)  -  Q)  (x!)  dx1  -  §0.  (-l)n(xj  +  Z-  —  j  sin  (C7) 

y  J  y  n=0  v  '  v  u0n  m=l  nut  a  '  v  1 ' 
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Figure  3.  Notation  for  stiffener  and  plate  forces 


Stiffener  temp.  =  0 

Plate  temp,  =  9  sin  —  sin 

a  a 

Coefficient  of  expansion  =  a 

(a)  Thermal-stress  problem  (b)  "Shear-lag"  problem. 

Figure  4.  -  Problems  considered  for  numerical  examples. 


Figure  5.  -  Variation  of  dimensionless  plate  and  stiffener  forces  for 
the  thermal-stress  problem,  (a)  X  =  tah/(rt2A)  =  0.0. 
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Figure  5.  ~  Continued,  (b)  \  =  4ah/ (jt2a)  =  O.lj  dashed  curves  represent 
approximate  solution  using  equations  (b8o). 


0-5 


Figure  5*  “  Concluded*  (c)  X  5  lj-ali/(:rr2A)  =  1*0 


to 


oo 


at  x/a  =  0,0 


Figure  6.  -  Variation  of  dimensionless  plate  and  stiffener  forces  for 
the  shear-lag  problem,  (a)  \  =  kah/(itsA)  -  0.0  . 


Continued,  (b)  \  =  hsii/ («SA)  =  0.1)  dashed  curves  represent 
approximate  solution  using  equations  (B8l)  to  (b84) . 


to  oo  at  x/a  =  0.0 


0.8 
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Figure  6.  -  Concluded. 

(c)  X  =  4ah/(naA) 

=  1.0  . 
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examples,  one  a  thermal- stress  problem,  the  other  examples,  one  a  thermal- stress  problem,  the  other 

a  "shear-lag"  problem.  a  "shear-lag"  problem. 


< 

to 

< 

Z 


to 

I 

a  s 

Xi  rO 
-M  O 

a  fe 


ID 

<D 

CO 

CJ  £ 

1 

Q 

°  Sf 

Eh 

CO  J 

5  a 

a  o 

as 

2 

(D  cti 

.8  g 

si 

*  & 
CD  _ 

§  Sf 

„  H 
CO  ■ 
CD  £ 
■—I  03 

a  <d 

a  m 


NASA  NASA 


CO  4J 

0  rt 

d 

x: 

O 

<D  rt 

o 

■°  5 


3*  *a  g 

Q  §  tD  H 

-  *«  0 
CO 


•T5  lo 

CM 

I 

Q 


•C  £ 
O  □ 

a  H 


rt  3  a  r-1 

Ami 

3  a  5  w  (S  s§ 


h  H  H  w 


s 

73  73  2  3 

0  0  cF  -£ 
73  XJ  Eh  « 

a  •£ 
a  5-1 
go 

■°  833 

45  q  o 
a°  ^2 
"S  ^  s 

r-H  -M  y 

P  O  5  D 
§b  0  3  > 
c  *P?  2  rt 


&  §  § 

X  °  o 

rt  Xj 

&  ® 
bp  ^  ? 
a 


0 


CD 

>>  > 
rt  a; 
73 


s  „ 
*&  ° 
Js  2  t? 

r.  Ctf  0 

H  r^-S 

&’C 

£  y 


SfSja 

CO  ^  o 

C  +J 


^  t;  cc  o 

fH  £ 

0  'S 

s  B 


O 

0 

i-t  73 

a 

o  rt 


73 

CO  CD 

c  s 


CO 

3  *  £ 

E  a  -a 

■m  id  P 

C0  J  X5 
_ ,  a  ,  i 
rt  ^ 
g  73 
O  <D 

1  g, 

I-2 

CD  a 

<D  ° 

■a  § 


if  g  -c 

<tt  3  rt  rt 

3  CD  O 

g  8  §  5s 

a  £  || 

■o  H  S  § 

CD  >  _ 

£  eg 

a  co  g  5 

CT1  CO  U 
ID  0)  M  o 
U  U  £»  ° 

+->  1—1  CQ 

•S  "  I 
«  8 
I  CD  +3  CD  a 
co  _,  ■jh  e 

M  »c  ^  . 

O  *}  o  £» 

>  CO  rQ 

2  g>.sl 

3.8  g  g 

CD  *  5  ^ 
X!  73  Pa  © 
-*->  CD  H 

a  e  o  rt 

~  s  a# 

.  a  e  a 

CO  ©  g  - 


-  a 
§■^1 
-y  73  73  2 


CO 

CO  a 
CD  CD 
C  X 


'  ^ 
§  S 

o  o 


a 

o  . 

CO  CO 
CD  73 
a  rt 
CD  O 


73  43  a 

£  04  ^  h 

§ 

^  I A  8 

.2t$  g-s 

CD  a 
CO  JL  73  rt 
•>h  C3 
®  o  rt  ® 

>1  •<-!  R 

d  w  o 
2  2  ^  £ 
5  5  <d  -a 
rt  aJ  c  k 

m  o  S3  h 


M— t  -f-H  ctf  4-> 

d  ID  _  U 

m  3J  §  m 
bp*°  0J  1 
e  S  co 

ia«s 

^S" 

2£  sf 

CD  H  o 

N  %u 

g  s  *  8 

£  I  -  2 

G  1-C  rv 

■“  a  -a 

S! 


s » 

f  I  ■§ 

<h  X  d 

o  a 


^  CD 

o  o 

CO  C 


§5  S>s 

-e  ^  a 


73  H 
CD 

fH 

*rH  • 

a  w 

o4  CO 
CD  CD 

u  u 

-M 
CD  CO 
JC 


CO 


CO 


CDXj  CO 
O  CO 
£h 

‘  O  bo  OJ 
2  X!  73  ^ 

St;  “5 

S°  s  2 

Cols, 


a  73 
O  CD 

is, 

a  a 

CD  ° 

I  § 

4h  O 

4H 

•s 

a 


a  £ 
o  ^ 
a  ^ 

5  Spa 

rt  o  ^ 

H  i-H 

cd 

-a  73  P>h 

^  a  o 
M  S 

,  a 

co  ^ 

CO  •S? 


«  O 

-  CO 

a 

*  rt 
CO  CD 

.s  a 

© 

CO  X5 


a  73 
CD 


s 


co  u]  -h  ,Q 


CO  X» 

=3 

ci 

CD  X5 

a  3 

>»  S 


|L  s 

ts-g-g^ 

2  *S  ^  H 
a  r 
a  a  h 
d  ,§  co  m  ^ 
(h  cd  a  c 

^3  B.8| 

®  a-H  >,  g 

£2  73  g  0 
d  CD  -o  N 

11  g 


CO  '  H  a  ra 

co  a  c  rt  0  .S 

CD  0  °  C  Qi 

a  ^  o  o  5*  g  73 

rt 

&  8 

2  CO 

*  a  (/j 

co  ■ 
co 
o 
a 
o 


U  ® 

CO  0 
bo  *a 


© 


<  o  ^ 
55  525  ^ 


a 

© 

CO 
© 

a  _ 

a  S  •“  ■o 

•2t5  S-g 

«  2  ■§  at 

»  o  d  « 

Is” 

“leg 

co  ©  Ph 
co  .2 .53  ij 
g  CDX3  co 

”  O  CO  Z3 

CO  ^3  0  73 
'  O  bJ3^ 
2  J3  73  U 

a  a  ^  -3 

ft  °  ^  a 

Colg, 


g73 

>*5 

CD 


rH  4J  O  rj 

'g  5  rt  rt 

°  3  CD  O 

{h  °  °  a 

I  "  g  s 
«J8  &S 

73  H  S 


CO  73 
cd  0 
.  X3 


iS 

o  +* 


rt  0 
a 
ft 


73 
CQ  0 

§  1 


>  co 

I  8 

Oigo 
^•H  3-ti 

+e  m  g 

CO  ,Q 

7j  -S  -{d 

a  a 

a  73  rt 

o  0  -£ 

to  -h  a  w 

rt  8  S.  § 


a  co 

O4  CO 

0  0  co  "5 

a  a  +*  ° 

-(->  M  ro 

S”  .g 

CO  0 

®  a 

a 

0  >> 
to 


0  0 

rt  £ 
CO 

a 

0 


|g 


5  g  '3 

CO  rD  a 


0 

CO  _, 

Ss 

S^l 
*f  sS 

XI  73  fa  0 

+-»  0  a 
el  o  rt 
m8»§ 

« .s  Ss 

!^SS 

au  .rt 
CO  X*  0 
10  -h  £!  <H 


"The  aeronautical  and  space  activities  of  the  United  States  shall  be 
conducted  so  as  to  contribute  ...  to  the  expansion  of  human  knowl¬ 
edge  of  phenomena  in  the  atmosphere  and  space.  The  Administration 
shall  provide  for  the  widest  practicable  and  appropriate  dissemination 
of  information  concerning  its  activities  and  the  results  thereof  ” 

— National  Aeronautics  and  Space  Act  of  1958 
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